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Abstract 

Vortex lattices are constructed in terms of linear combinations of solutions for 
Scrodinger equation with a constant potential. The vortex lattices are mapped on 
the spaces with two-dimensional rotationally symmetric potentials by using confor- 
mal mappings and the differences of the mapped vortex-patterns are examined. The 
existence of vortex dipole and quadrupole is also pointed out. 
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Vortices play an interesting roles in various aspects of present-day physics such as vor- 
tex matters (vortex lattices) in condensed matters [1,2], quantum Hall effects [3-7], various 
vortex patterns of non-neutral plasma [8-11] and Bose-Einstein gases [12-16]. The vortices 
are well known objects in hydrodynamics and has been investigated in many aspects [17- 
21]. In quantum mechanics hydrodynamical approach was vigorously investigated in the 
early stage of the development of quantum mechanics [22-29] . The fundamental properties 
of vortices in quantum mechanics were extensively examined by Hirschfelder and others 
[30-34] and the motions of vortex lines were also studied [35,36]. Recently an analysis 
of vortices has been carried out on the basis of non-linear Schrodinger equations [37,38]. 
The vortex dynamics will hold an important position in quantum mechanical problems. 
In vortex dynamics, however, the construction of various vortex patterns from fundamen- 
tal dynamics is still not very clear. Recently Kobayashi and Shimbori has proposed a 
way to investigate vortex patterns from the special solutions of Schrodinger equations 
[39]. By using the conformal mappings the article has shown that the special solutions 
with zero energy eigenvalue in the two dimensional parabolic potential barriers (2D PPB) 
expressed hy V — — 777,7^ (x^ + y^)/2 [40] can be extended to all rotationally symmetric 
potentials and the solutions are infinitely degenerate as same as those in the 2D PPB [39]. 
It has also been shown that the infinite degeneracy can be an origin of variety of vortex 
patterns which appear at nodal points of wave functions and this idea can be applied in 
three dimensions. In the article a few examples of simple vortex patterns are presented. 
In this letter some patterns of vortex lattices that are constructed from simple linear 
combinations of the infinitely degenerate states will be studied. 

Let us start fom the arguments in two dimensions [39]. In two dimensions the eigen- 
value problems with the energy eigenvalue S are explicitly written as 

[-^ A +K(p)] ^{x, y)=E ^(x, y), (1) 
where A= d"^ jdx^ -|- d"^ /dy'^, the rotationally symmetric potentials are generally given by 

Va{p) = -a'gy^^-'\ 

with p = a/x^ + y^, a G M (a 7^ 0), and m and Qa are, respectively, the mass of the particle 
and the coupling constant. Note here that the eigenvalues S should generally be taken as 
complex numbers that are allowed only in conjugate spaces of Gel'fand triplets [41]. Note 
also that T4 represents repulsive potentials for {ga > 0, a > 1) and {ga < 0, a < 1) and 
attractive potentials for {ga > 0, a < 1) and {ga < 0, a > 1). 
Let us consider the conformal mappings 

= z", with z ^ x + iy. (2) 

We use the notations and Va defined by Ca = 'i*a + iVa that are written as 

Ua — cos a(p, Va — siu a(f, (3) 



1 



where (/? = arctan(y/a:). In the {ua,Va) plane the equations (1) are written down as 

where Aa= + d'^ ldv\. We can rewrite the equations as 

[-^ A„ -g^\ ^K, v:) = a-'£ pf -"^/'^^K, Va). (5) 

It is surprising that the equations become same for all a 7^ when the energy eigenvalue 
S have zero value. That is to say, for £^ = the equations have the same form as that for 
the free particle with the constant potentials ga as 

[-7r~ ^« ~9a] l/^{Ua,Va) = 0. (6) 

zm 

It shoTild be noticed that in the case of a = 1 where the original potential is a constant gi 
the energy does not need to be zero but can take arbitrary real numbers, because the right- 
hand side of (5) has no p dependence. In the a = 1 case, therefore, we should take gi + £ 
instead of ga- Here let us briefly comment on the conformal mappings Ca — z"'- We see that 
the transformation maps the part of the {x, y) plane described by < p < oo, < </? < tt/o 
on the upper half-plane of the Vg) plane for a > 1 and the lower half-plane for a < —1. 
Note here that the maps on the part of the {ua, Va) plane with the angle (fa = (f + o. can 
be carried out by using the conformal mappings 

Ca{a) = z"e-. (7) 

In the maps the variables are given by 

Ua{a) = UaCosa — VaSma, Va{a) = VaCosa + UaSma. (8) 

The relations Ua{0) — Ua and Va{0) — Va are obvious. 

It is trivial that the equations for all a have the particular solutions 

i^^iua) = Nae^'^^^^ (9) 

^o^K) = iVae^''""% for^„>0 (10) 

and 

4>^{ua) = M^e^'^""" (11) 

0±(^;„)=M„e±'=«''", for^„<0 (12) 

where ka = ^J'lm\ga\lh,^ and A^^ and are in general complex numbers. General solu- 
tions should be written by the linear combinations of (9) and (10) for gi^ > and those 
of (11) and (12) for g^ < 0. Examples for the 2D PPB with ga > are presented in ref. 
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[40] . In the following investigations we shall concentrate our attention on the solutions of 

(9) and (10) that are expressed in terms of plane waves. Hereafter we use the notations 
Ua{a) and Va{a) with — tt < o; < tt, which are used in the conformal mappings of (7). We 
see that 

are the solutions of (6). (For details, see ref. [39].) 

We have to comment on the degeneracy of the solutions. The origin of the infinite 
degeneracy can easily understand in the case of the 2D PPB [40]. It is known that 
energy eigenvalues of ID PPB are given by pure imaginary values ^i{n + 1/2)^7 with 
n = 0,1, 2,--- [42-47]. From this reult we see that the energy eigenvalues of the 2D 
PPB, which are composed of the sum of the ID-PPB eigenvalues with the oposite signs 
such as i{nx — ny)hj with Ux and Uy = 0, 1, 2, • • • , include the zero energy for n^, = Uy. 
It is apparent that all the states with the energy eigenvalues i{nx — ny)hrf are infinitely 
degenerate [40]. This means that all the rotationally symmetric potentials have the same 
degeneracy for the zero energy states. By putting the wave function f^{ua;Va)4'^{ua) 
into (6) where f^{ua,Va) is a polynomial function of Ua and Va, we obtain the equation 

[^a^'^^K^]f^{Ua,Va)=Q. (13) 
dUa 

A few examples of the functions / are given by [40] 

1, 

A{Aklvl + l±AikaUa). (14) 

Wc can obtain the general forms of the polynomials in the 2D PPB, which are generally 
written by the multiple of the polynomials of degree n, H^[\/2k2x), such that 

f^{u2;v2) = H^{V2hx) ■ H^{^2y). (15) 

where x and y in the right-hand side should be considered as the functions of U2 and V2 
[40]. Since the form of the equations (6) is the same for all a, the solutions can be written 
by the same polynomial functions that are given in (15) for the PPB. That is to say, we 
can obtain the polynomials for arbitrary a by replacing U2 and V2 with Ua and Va in (15). 
Note that the polynomials H^{^) with ^ = \/m/yJhx are defined by the solutions for the 
eigenstates with S = ^i{n + l/2)hj in ID PPB of the type V{x) — —7717^x^/2 and they 
are written in terms of the Hermite polynomials i?n(0 

(For details, see ref. [45,46].) The states expressed by these wave functions belong to the 
conjugate spaces of Gel'fand triplets of which nuclear space is given by Schwarz space. 
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Actually we easily see that the wave functions cannot be normalized in terms of Dirac's 
delta functions except the lowest polynomial solutions [39]. 

The extension to three dimensions can easily be carried out in the cases with potentials 
that are separable into the {x, y) plane and the z direction such that 

V{x,y,z)^Va{p) + V{z). 

When the energy eigenvalues of the z direction is given by Ez, we obtain the same equation 
as (6) for £ — Ez = (). (For the a = 1 case gi + £ — Ez > Q should be taken.) If we 
take the free motion with the momentum Pz for the z direction, E^ = pl/2m should be 
taken. It is important that the total energy £ is in general not equal to zero in the three 
dimensions. Note that wave functions for the separable potentials are written by the 
product such as ip{x,y,z) — ip{x,y)ip{z). Hereafter we shall not explicitly write ip{z) in 
the wave functions. 

Let us study vortices that appear in the hnear combinations composed of the solutions 
with the polynomials of (15). Before going into the details we briefly describe vortices in 
quantum mechanical hydrodynamics. The probability density p{t, x, y) and the probabil- 
ity current j{t,x,y) of a wave function ip{t,x,y) in non-relativistic quantum mechanics 
are defined by 

p(t,x,y) = \ip{t,x,y)f , (17) 
j{t, X, y) = Re [iP{t, X, y)* {-ifiV) iP{t, x, y)] /m. (18) 

They satisfy the equation of continuity dp/dt + V • j = 0. Following the analogue of the 
hydrodynamical approach [17-21], the fluid can be represented by the density p and the 
fluid velocity v. They satisfy Euler's equation of continuity 

^ + ipv) = 0. (19) 

Comparing this equation with the continuity equation, we are thus led to the following 
definition for the quantum velocity of the state ijj{t, x, y), 

V ^ Mr^yL^, (20) 

\ij{t,x,y)\ 

in which j{t, x, y) is given by (17). Notice that p and j in the present cases do not depend 
on time t. Now it is obvious that vortices appear at the zero points of the density, that 
is, the nodal points of the wave function, where the current j does not vanish [39] . We 
should here remember that the solutions of (6) degenerate infinitely. This fact indicates 
that we can construct wave functions having the nodal points at arbitrary positions in 
terms of linear combinations of the infinitely degenerate solutions [39,40]. 

The strength of vortex is characterized by the circulation F that is represented by the 
integral round a closed contour C encircling the vortex such that 

r = j v-ds (21) 
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and it is quantized as 



r = 2T:lh/m, (22) 

where the circulation number I is an integer [31,33,36]. It should be stressd that we 
can perform the investigation of vortices for all the cases with a 7^ in the fa) 
plane, because fundamental properties of vortices such as the numbers of vortices in the 
original plane and the mapped plane and the strengths of vortices do not change by the 
conformal mappings. Vortex patterns in the {x, y) plane can be obtained by the inverse 
transformations of the conformal mappings. Let us here show that vortex lines and vortex 
lattices can be constructed from simple linear combinations of the low lying polynomial 
solutions. And also the mapped patterns of those lines and lattices are investigated by 
the conformal mappings for the a = 2 (PPB case; Va oc p^) and for the a = 1/2 (Coulomb 
type; oc p~^). In the following discussions the suffices a of Va and ka are omitted. 

(I) Vortex lines 

Let us consider the linear combination of two degree 1 solutions such that 

^{u,v)^v(^^'' -ue-'^\ (23) 

where the complex constant corresponding to the overall factor of the wave function is 
ignored, because the wave function belongs to the conjugate space of Gel'fand triplet and 
it is not normalizable. This means that the wave function represents a stationary flow 
[40] . The nodal points of the probabihty density 

|\l'(M,f)| =u +v — 2uv COS k{u + V) 

appear at points satisfying the conditions 

u = ±v, cos k{u + v) = ±1. (24) 

We have the nodal points at 

u = v = mr/k, for n integers. (25) 

In the {u, v) plane the positions of vortices can be on a line ofu = v. After some elementary 
but tedious calculations we see that the circulation numbers of vortex strengths are given 
by / = —1 for n =positive integers and / = 1 for n =negative ones. Note that the 
origin a,t u — v — has no vortex. We can directly see the result by showing the fact 
that the strength of vortex F becomes zero for the closed circle around the origin. We 
can interpret this result as follows; at the origin there exist a pair of vortices having the 
opposite circulation numbers, that is, they, respectively, belong to the vortex line with 
I — —1 and that with 1 = 1. We may say that it is a vortex dipole. 

For the case of a = 1 (constant potential) we can take as u = x and v = y. 
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In the case oi a — 2 (PPB) we have u — — and v — 2xy. The relations for the 
nodal points are written down as 



1 , 1 I nn . ... 

y — — -X, y — ^-^j-;—^ — for n = positive integers 




1 , 1 / n TT _ 
y — ;= X, y — ±-A / — , for n — negative integers. 

(26) 

We see that a vortex quadrupole composed of two vortex dipole appears at the origin. 

In the case of a = 1/2 (Coulomb type), by using the relations — v^ — x and 2uv = y, 
we obtain the conditions for the nodal points as follows; 

X — 0, y = 2 , for n — non zero integers. (27) 

Note that the origin is a singular point, where the source of the potential exists. 

Figures for a = 1,2 and 1/2 are presented in figs.l, 2 and 3, which, respectively, 
represent the vortex pattern for the constant potential, that for the PPB, and that for 
the Coulomb type one. Note here that the differences of the potentials clearly appear 
not only in the vortex patterns but the properties of the singularities at the origin as 
well. In hydrodynamics it is known that vortex lines for the constant potential (a = 1) 
are stationary but generally unstable for perturbations [17-21]. Note that the parallel 
vortex lines are constructed from the linear combinations of the lowest and the degree 1 
polynomials 



(II) Vortex lattices 

Let us consider the linear combination of a stationary wave and a plane wave such 
that 

'^(u,v) ^cosku-e-'''\ (28) 
The nodal points of the probability density 

\^{u,v)\^ — 1 + cos'^ ku — 2coskucoskv 
appear at positions satisfying 

u — rmr/k, v — mr/k, (29) 

where both of m and n must be even or odd, that is, (—1)" = (—1)" must be ful- 
filled. These conditions produce a vortex lattice presented in fig.4, which was suggested 
in ref.[48]. 
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In the cases of the PPB (a = 2) and the Coulomb type {a — 1/2) vortices appear at 
the cross points of the following two functions; 

-y^ ^ mn/k, xy = m:/2k, for the PPB, 

+ = (m^ + n^fir'^/k'^, y = 2mr^7^7A;^ for the Coulomb type. (30) 

In these case we obtain the circulation number I — —1 for the all vortices. Figures for 
a — 2 and 1/2 are given in figs. 5 and 6, respectively. 
In these arguments we see the following points: 

(1) The construction of vortex lattices in experiments seems to be not very difficult. In 
fact the vortex lattice of (II) can be produced from a stationary wave and a plane wave 
perpendicular to the stationary wave. 

(2) The differences of potentials can be clearly seen from the vortex patterns. Especially 
the distances between two neighbouring vortices are a good object to identify the type of 
the potentials. That is to say, the vortices appear in an equal distance Ti/k in the case of 
the constant potential (a = 1), whereas the distances become smaller in the regions far 
from the origin for a > 1 and larger for a < 1 in comparison with those near the origin. 

(3) The property of the singularity at the origin is also a good object to identify the 
potentials. 

Though it is at this moment difficult to categorize the present experimental vortex patterns 
[8-16], we shall be able to understand fundamental dynamics of vortex phenomena from 
vortex patterns. 

It is also noticed that the present results can be applicable not only to quantum 
phenomena but also those in classical fluids by changing the parameters m, h and ga in 
the original equation. 

Up to now we have not discussed on the stability of the vortex lattices. In order 
to investigate the time development of the patterns we have to take account of the fact 
that the solutions used here belong to the conjugate spaces of Gel'fand triplets. In the 
spaces the eigenstates generally have complex energy eigenvalues such as resonances and 
the eigenvalues arc expressed by pairs of complex conjugates such that S = e ± 27, where 
e, 7 G M [41]. It is known that the + and - signs of the imaginary part in the eigenvalues, 
respectively, represent resonance-formation and resonance-decay processes. We see that 
this pairing property is the origin of the infinite degeneracy of the solutions and the infinite 
degeneracy stems from the balance between the resonance-formation and resonance-decay 
processes. This fact seems to indicate that many vortex systems are possibly unstable for 
perturbations. Actually the existence of vortex lattices has already been pointed out and 
it has also been noticed that those systems will decay from their edges, where the balance 
between the formation and decay processes is broken [48]. In such processes we have to 
discuss on dynamics in many vortex systems and statistical mechanics should be extended 
into the conjugate spaces of Gel'fand triplets, where the freedom of the imaginary energy 
eigenvalues must be introduced [48-50]. In the theory the infinite degeneracy plays another 
important role, that is, it becomes the origin of new entropy. The new entropy brings 
essentially new processes in thermal non-equilibrium [49,50]. At present we still have 
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many fundamental questions in the investigation of dynamics in the conjugate spaces of 
Gel'fand triplets. The study of vortex dynamics will brings us a new prospect and open 
a new dynamics for the systems essentially described by the states in Gel'fand triplets. 
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Fig. 1: Positions of vortices for n = ±1, ±2, ±3 in the constant potential (a = 1), which 
are denoted by •, and o stands for the vortex dipole at the origin. 




Fig. 2: Positions of vortices for n — ±1, ±2, ±3 in the PPB {a — 2), which are denoted 
by •, and o stands for the vortex quadrupole at the origin. 
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Fig. 3: Positions of vortices for n = 1, 2, 3 in the Coulomb type potential (a = 1/2), 
which are denoted by •, and o stands for the source at the origin. 
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Fig. 4: Positions of vortices for m, n = 0, ±1, ±2, ±3, ±4 in a constant potential 
(a = 1), which are denoted by • and the distance between the neighbouring hnes are 
taken by ir/k. 
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